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吴 差 环境 中 参数 识 辨 前 测量 信息 的 灶 描 述 
李 静 !， 丁 海洋 !， 任 学 部 


(1. 国防 科技 大 学 信息 通信 学 院 , 西安 710106; 2. 空军 工程 大 学 , 西安 710043) 


摘 要 : 在 逆 问 题 的 参数 识 准 中， 测量 信 息 包含 的 信息 量 直 接 影响 参数 的 重 构 精 度 ， 如 何 度 量 测量 信息 中 的 信息 量 ， 
对 选取 参数 识 辨 所 需 的 测量 点 个 数 具有 重要 的 指导 人 作用。 根据 贝 叶 斯 方法 ， 将 体现 先 验 信息 的 先 验 概率 和 反映 测量 信 
息 的 似 然 概率 融合 得 到 待 识 办 参数 样本 出 现 的 概率 ， 结 合 信 息 米 和 最 大 业 计 算 评价 因子 ， 将 评价 因子 用 于 描述 参数 识 
状 前 带 有 误差 的 测量 信息 所 包含 的 信息 量 。 计 算 实 例 表 明 ， 该 计算 方法 可 以 有 效 地 描述 误差 环境 中 的 测量 信息 ， 对 实 
际 应 用 中 测量 信息 的 选取 具有 重要 的 指导 意义 。 
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Entropy description of measured information with error before parameter identification 


Li Jing!, Ding Haiyangl, Ren Xueyao’ 
ls (1 College of Information & Communication, National University of Defense Technology, Xi’?an 710106, China; 2. Air Force 
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Abstract: In the parameter identification of inverse problems, it was well-known that the information quantity contained in 


measurement information influences the reconstruction precision of parameter directly. How to describe the information quantity 


in measurement information plays an important role in the selection of the number of measurement points. According to the 


Bayesian method, the probability of the parameter sample to be identified is calculated by combining the priori probability 


I reflecting prior information and the likelihood probability reflecting measurement information. The evaluation factors are 
calculated by the maximum entropy and information entropy of the parameter. This paper introduced the information entropy to 
describe the information quantity of measurement information with error before parameter identification. Numerical tests 
showed that this computing method proposed could describe the measured information with error effectively and determine the 
measured information in practice. 
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0 ”引言 效 的 数学 方法 一 一 测量 信息 的 炳 描述 由。 
与 过 去 文献 中 灶 的 用 法 E9 相 比 , 测量 信息 的 灶 描 述 将 箭 用 
参数 识 辩 是 逆 问 题 中 的 一 类 重要 问题 ， 在 逆 问 题 的 工程 应 ”在 参数 识 辩 前 测量 信息 的 描述 中 内 。 本 文 针 对 参数 识 辨 前 大 量 
用 中 ,参数 识 辨 前 会 提供 大 量 的 测 和 带 有 误差 的 测量 信息 ， 将 信息 炉 用 于 描述 测量 信息 所 包含 的 信 
为 了 提高 参数 识 辨 的 精确 度 ,所 用 的 测量 信息 越 多 越 好 ;但 是 ， 息 量 , 结合 贝 叶 斯 方法 "1, 采用 待 识 辨 参数 分 布 的 后 验 概率 定 
在 误差 环境 中 ， 这 些 测 量 信 息 往往 含有 一 定 的 测量 误差 ， 其 数 。” 义 参 数 样 本 出 现 的 概率 ， 计 算出 误差 环境 中 不 同 数量 的 测量 


量 的 增多 不 仅 会 加 大 计算 成 本 ， 更 重要 的 是 会 导致 重 构 计算 中 息 对 应 的 评价 因子 ， 通 过 比较 评价 因子 和 重 构 误差 ， 说 明 评价 
误差 的 积累 。 所 以 ， 在 逆 问 题 的 计算 过 获得 测量 信息 。 因子 可 以 用 来 判断 参数 识 辨 前 实际 所 需 的 测量 信息 ， 解 决 逆 问 
和 后， 通常 需要 先 对 大 量 数据 进行 第 选 ， 即 进行 参数 识 辨 前 大 题 应 用 中 的 首要 问题 。 不 失 一 般 性 ， 本 文 在 最 后 引入 一 实例 对 
量 数据 的 分 析 和 评 作 的 压力 和 重 构 误 差 。 误差 环境 中 测量 信息 的 描述 方法 进行 说 明 。 计 算 结果 表明 ， 
在 参数 识 辩 问题 的 研究 中 ， 已 有 的 大 量 文献 主要 针对 重 构 方法 ”以 根据 评价 因子 较 好 地 进行 逆 问 题 中 带 误差 的 测量 信息 的 描 
进行 了 讨论 0 53, 而 忽略 了 重 构 前 必须 解决 的 重要 问题 一 一 测量 
言 息 的 选取 。 目 前 ， 针 对 这 一 问题 本 文 作者 已 经 提出 了 一 种 有 
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录用 入 李 静 ,等 : a¥ a a YE 
1 ”评价 因子 的 定义 以 及 测量 信息 的 描述 HPO Dn ron (7) 
1.1 最 大 灼 确定 为 了 描述 似 然 概率 的 本 质 ， 

在 信息 论 中 ,信息 炉 就 是 随机 事件 不 确定 性 大 小 的 度量 09。 ne [4] (8) 
若 道 问题 的 数学 模型 用 算 子 方程 表示 为 :Af = y，f 是 待 识 辩 而 误差 六 《>1 时 ) 可 由 算 子 方程 左右 两 端 差 值 的 范 数 给 
参数 ，y 是 测量 项 ， 它 关于 f 的 有 限 次 测量 值 是 yy ,yy,……,》,， 出 [1; 

4 是 将 /映射 到 测量 值 7,y》,…, ,的 正 向 算 子 。 如果 随 机 事件 _]4- 川 而 
8 具有 多 种 可 能 解 fi, 万 …, ,每 个 解 出 现 的 概率 为 B,PB…,PP， m-—2 
则 该 随机 事件 对 应 的 信息 粹 为 : 综合 式 (8) 与 (9) ， 代 入 式 (7) 

H(f) Fana- dens (1) H.,(f) De ine (10) 

但 是 ， 如 何 量化 测量 信息 对 的 影响 呢 ? 根据 贝 叶 斯 方法 ， 比 式 即 为 逆 问 题 中 已 知 闫 个 测量 数据 后 参数 了 信息 粒 的 最 
本 文采 用 随机 事件 的 后 验 概率 来 定义 已 中 这样, 就 可 以 得 到 终 表达 式 。 
参数 了 在 已 知 不 同 测量 信息 后 的 信息 粹 。 通常 对 于 实际 问题 中 测量 信息 的 选取 ， 必 须要 对 测量 信息 

根据 贝 叶 斯 假设 中, 当 无 任何 先 验 信息 时 , 参数 1 可 能 的 解 ”给 出 统一 的 数学 方法 进行 描述 内。 定义 ， 

;了 太 …,Jy 对 应 的 先 验 概率 元 (f) (i=42…,N ) 应 在 /的 取 7=m| -HD) ei 
值 范围 D 内 是 “均匀 ”分 布 的 ， 所 以 z( 了 ) 可 写成 Han (f) 
为 已 知 m 个 测量 数据 后 的 评价 因子 。 它 可 以 在 待 识 辨 参数 
TE 入 当 / eD 12..N 四 值 未 知 的 情况 下 ， ee 
0 feD 评价 ， 并 可 以 对 参数 识 辨 中 实际 所 需 的 测量 信息 进行 分 析 。 
一 般 情况 下 ， 参 数 f 的 样本 空间 入 大 ， 不 能 对 所 有 样本 进 
其 中 :V 为 参数 的 样本 数 。 本文 令 吾 ,(7) 为 1 对 应 的 最 大 粹 (无 行 逐 一 计算 ， 所 以 采用 马尔 可 夫 链 蒙特 卡 罗 法 一 一 Gibbs 抽样 
任何 先 验 信息 时 的 炉 ) 外 ， 则 五 (了) 为 器 对 总 样本 抽样 ， 产 生 一 组 服从 后 验 概率 五 ( 门 六 ,yo) 分 布 的 
y Ny] 1 样本 。 下 面 用 一 个 具体 实例 进行 分 析 说 明 上 述 方法 的 有 效 性 。 
Hr (f) =D n= TN @) 
N 2 ”实例 分 析 
1.2 评价 因子 及 测量 信息 的 描述 方法 [4] 利用 一 维 的 参数 识 辨 模型 所 来 说 明 利用 上 述 方法 进行 测量 
贝 叶 斯 方法 的 基本 思想 是 将 先 验 概率 和 似 然 概率 融合 ， 信息 的 选取 。 该 模型 为 
到 问题 的 后 验 概率 。 根 据 贝 叶 斯 方法 的 原理 ， 有 贝 叶 斯 公式 : ey 
+ f(x):y(x)=0 xe[-1,1] (12) 
zi(f yy ) = RO yf) y(2) =1.71653 | 
Da pa |f) ® 
BS 其 中 :Fo 为 待 识 状 参数 ，y(») 为 测量 项 。 当 .Fo 的 解析 式 已 
式 (4) 右 端的 分 母 是 一 个 常数 , 于 是 元 (站 my) 正 比 于 分 。 知 : 
于 : f(x)=cos(x) —sin’ (xX) xe[—1,1] 
AGA ADLAG (5) (13) 
其 中 : f 是 待 识 辨 参数 ，z(f) 是 了 的 先 验 概率 ， 体 现 了 其 先 验 时 , 结合 式 (12) 可 以 得 到 yo 的 真实 分 布 , 如 图 1 所 示 。 
信息 ; yy,…,y 是 在 误差 环境 中 获得 的 测量 信息 ; a nd 
(94,…,y| 了 ) 是 包含 了 测量 信息 的 似 然 概率 ，(f 3%.…,y,) 是 训 ~ 
/的 后 验 概率 ， 用 它 来 表示 式 〔1) 中 的 只 ， 它 综合 了 参数 /的 本 \ 
先 验 信息 和 yy,…,y, 提供 的 测量 信息 ， 反 映 的 是 参数 了 的 综合 22 / * 
信息 。 该 后 验 概率 定义 在 整个 解 空间 上 , 表示 问题 的 “完全 ” 解 ， i / 
就 是 由 所 有 可 能 出 现 的 解 组 成 的 解 集中 任何 一 个 解 出 现 的 概 / \ 
率 。 证 四 


| 


将 式 (2) 与 (5) 结合 ， 式 (5) 可 简化 为 1 y(x) 的 真实 解 
Tf ly sy,) TVs | 亡 
\ 彤 为 


此 时 ， 参 数 / 对 应 的 信息 


(6) 


本 文 将 YD 沿 x 轴 在 [-1,1] 内 按 1/8 间隔 均匀 离散 ， 共 得 到 
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17 个 测量 数据 ， 给 这 些 测量 信 ， 
15% 的 相对 误差 后 ， 将 其 
“CX)a 


和 和 土 


T= 


> 


下 
将 ) 


息 分 别 加 入 +3%、 土 5%、 土 10% 


为 了 验证 前 


f(x) 的 重 构 误 差 : 


信 
低 ， 则 重 构 所 | 


Ax=1/8， 则 离散 


再 将 重 构 误差 与 计 
能 够 代替 重 构 误 差 较 好 地 
息 。 评 价 因子 越 大 ， 


54LO-/ OW |x100% 
f(D 


作为 不 同 误差 环境 中 的 测量 信息 


掉 关 于 逆 问 题 参数 识 辩 前 测量 信息 的 描述 方法 
四 本 文 假设 f(x) 的 解析 式 未 知 。 已 知 y(x) 的 测量 信息 y*(x)， 
限 差分 法 求 得 的 .Fo 的 解 f°(x) 与 其 解析 解 


FE 价 因子 进行 对 比 ， 从 而 说 明 用 评价 因子 
述 不 同 误差 环境 中 重 构 所 需 的 测量 
i ne 


与 ?00 相同 ， 


值 范围 六 = [5.0， 再 将 /Go 在 [-1/5 汕 内 按 Af =2/15 的 间隔 等 
间隔 离散 ， 的 个 数 1 = 10 。 则 两 次 离散 后 AD 的 总 


样 


限 差分 对 全 只 展 
dx 


其 中 ,加 是 测量 数 

对 于 测量 信息 y 
求 得 M 个 + 值 ， 
数据 对 应 的 评价 因子 ， 结 


J (Co 重 构 误差 的 平均 。 
图 2 和 3 中 测量 点 的 个 数 是 按照 图 1 中 y (2) 的 离散 分 量 
2 可 以 看 出 ， 在 没有 测量 误差 的 情 
也 反映 重 构 误差 的 变化 规 得 
[ 重 构 误差 都 越 小 。 

不 难看 出 ， 当 加 入 误差 之 后 ， 已 知 
重 构 误差 和 评价 因子 并 不 一 定 越 小 
构 误差 和 评价 因子 并 不 是 最 小 )。 这 主要 
息 进行 重 构 时 ， 在 测量 信息 增多 


从 左 至 右 依次 递增 。 
况 下 ， 评 价 因子 能 够 较 好 二 
点 数 越 多 ， 评 价 

观察 图 3 的 (a) 
越 多 ， 其 对 应 也 
17 个 测量 点 时 的 习 
是 因为 ， 利 用 带 误 差 的 测量 


的 测量 信 


( 


无 


的 


还 显示 ， 本 


的 测量 


则 离散 分 二 
本 数 W= 上 =107 。 
本 例 中 ， 随 机 抽取 样 


息 不 足 或 不 够 准确 。 
在 [-1,1] 内 均匀 离散 ， 离 散 间隔 
量 的 个 数 k=17 。 根 据 先 验 信息 ， 了 (%) 的 


以 
A 


去 


本 数 M =2x10;， 根 据 式 (3) ，f(%9) 
对 应 的 最 大 炉 矿 ,(f)=InM =In2x10” =12.2061。 然后 运用 有 


将 式 (12) 代入 式 (9) ， 此 时 一 即 为 


-2700)+y(G0 -2A0]/ (2A) + fe)- ye) 


将 TT 


m-2 (15) 
i=1,,M; j=,m 


别 将 f() 的 MM 个 样本 代入 式 (15) 
人 入 式 〈10) (11)， 便 可 得 到 m 个 测量 
图 2 和 3 所 示 。 由 于 已 知 测量 信 
息 是 离散 的 ， 图 2 和 3 中 的 重 构 误差 是 在 yc 相应 测量 点 上 


证 


中 


， 并 且 测 


比较 ， 有 参数 


关 


重 构 误 关 


44) 


息 的 个 数 是 5。 与 17 个 
点 进行 重 构 可 以 节约 一 定 


0.08 


oor] 


ooe 


甚至 更 多 的 测量 点 相 比 ， 利 月 
E 的 计算 成 本 。 


0 一 评价 因子 
[4 一重 构 误差 ] 


a S A 
0—0=—8—A<0 


测量 信和 点 的 个 数 


2 无 测量 误差 时 ,不 同 数量 的 测量 信息 对 应 的 评价 因 


10 12 14 16 
测量 信息 的 个 数 (a) 


T T T 
10 12 14 16 
测量 信息 的 个 数 (b) 


此 ， 在 逆 问 题 包 


3 不 同 的 误差 环境 中 ， 不 同 数量 的 测量 信息 对 应 的 评 
误差 


同时 ， 其 携带 的 误差 有 可 能 时 致 重 构 误 差 的 积累 。 同 时 图 3 


构 误差 和 评价 因子 均 会 随 着 相对 误差 的 增加 而 增 大 ， 


但 是 ， 当 相对 误差 不 断 变 大 时 , 图 (b) 中 重 构 误差 的 变化 越 来 


越剧 烈 ， 


以 得 出 这 样 的 结论 


反 


过 


而 图 (a) 中 鹿 


图 3 可 以 得 到 这 样 世 
因子 最 小 和 重 构 误 差 最 低 的 时 候 ， 


价 因子 的 变化 却 是 越 来 越 平缓 。 由 此 可 

: 当 相 对 误差 过 大 时 ， 评 价 因 子 无 法 很 好 地 

映 重 构 误差 地 变化 规律 ， 不 能 用 来 进行 测量 信息 的 描述 。 通 

论 : 测量 信息 的 个 数 为 5 时， 是 评价 
天 


此 ， 该 重 构 需要 的 测量 信 


对 误差 环境 中 的 测量 信 | 


识 辨 ， 又 可 以 贡生 
结束 语 


在 逆 问 题 的 工程 应 
题 。 本 文 针对 实际 应 用 中 的 这 一 问题 ， 


[ 程 应 用 中 ， 参 数 识 辩 前 利用 评价 
息 进行 描述 ， 既 可 以 较 准确 
大 量 的 计算 成 本 。 


寺中 ， 测 量 信 


评价 因子 代替 重 构 误差 对 误差 环境 


了 描述 。 同 时 , 将 信息 炉 引 入 到 逆 问 题 
We ee 


介 因 子 和 重 


筷 的 选取 是 一 


参数 识 辩 前 


,应 月 


描述 为 先 验 概率 ， 测 量 信 


然 概率 ， 二 者 


， 等 : 误差 环境 中 参数 识 准 前 测量 信息 的 精 描 术 


重 构 误 差 


个 必须 解决 
首次 提出 如 何 
上 量 信息 


里 论 ， 
息 表示 为 似 


到 已 知 不 同 数 量 的 测量 


参数 的 最 大 糯 和 信息 
因子 。 在 本 文 的 最 后 ， 


用 信 


息 粮 表示 的 评价 


a i de 
乱 时 参数 的 信息 业 。 
gi 
体 


用 待 识 辨 
F 价 
的 实例 分 析 表 明 : 相对 误差 不 


因子 能 够 较 好 地 反映 重 构 误差 
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的 变化 规律 ， 可 以 用 来 完成 逆 问 题 中 测量 信息 的 描述 。 此 方法 
不 仅 能 满足 重 构 误差 低 的 要 求 ， 还 可 以 大 大 降低 计算 量 。 该 项 
研究 对 逆 问 题 的 工程 应 用 具有 重要 的 指导 意义 。 
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